Fractal Weyl law for quantum fractal eigenstates 
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The properties of the resonant Gamow states are studied numerically in the semiclassical limit for 
the quantum Chirikov standard map with absorption. It is shown that the number of such states is 
described by the fractal Weyl law and their Husimi distributions closely follow the strange repeller 
set formed by classical orbits nonescaping in future times. For large matrices the distribution of 
escape rates converges to a fixed shape profile characterized by a spectral gap related to the classical 
escape rate. 
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The Weyl law [1] gives a fundamental link between the 
properties of quantum eigenstates in closed Hamiltonian 
systems, the Planck constant h and the classical phase 
space volume. The number of states in this case is de- 
termined by the system dimension d and the situation is 
now well understood both on mathematical and physical 
grounds [2, 3]. Surprisingly, only recently it has been 
realized that the case of nonunitary operators describing 
open systems in the semiclassical limit has a number of 
new interesting properties and the concept of the fractal 
Weyl law has been introduced to describe the dependence 
of number of resonant Gamow eigenstates on fi [4, 5]. 
The Gamow eigenstates find important applications for 
decay of radioactive nuclei [6], quantum chemistry re- 
actions [7], chaotic scattering [8] and microlasers with 
chaotic resonators [9, 10, 11]. Thus the understanding of 
their properties in the semiclassical limit represents an 
important challenge. 

According to the fractal Weyl law [4, 5] the number of 
Gamow eigenstates 7V 7 , which have escape rates 7 in a 
finite band width < 7 < 75, scales as 

-(d-i) 



N 7 oc h 



(i) 



where d is a fractal dimension of a classical strange re- 
peller formed by classical orbits nonescaping in future 
(or past) times. By numerical simulations it has been 
shown that the law (1) works for a 3-disk system [5] and 
quantum chaos maps with absorption [12, 13] at spe- 
cific values of d. Recent mathematical results for open 
quantum maps are presented in [14]. The law (1) for 
open systems with a fractal dimension d < 2 leads to 
a striking consequence: only a relatively small fraction 
of eigenstates /i ~ N 1 /N ex TS 2 ~ d>} have finite values of 
7 while almost all eigenstates of matrix operator of size 
N oc 1/h have infinitely large 7. The later states are as- 
sociated [12] with classical orbits which escapes from the 
system after the Ehrenfest time [15]. The former states 
with finite 7 are related to the classical fractal repeller 
and have been named quantum fractal eigenstates due to 
a fractal structure of their Husimi distributions closely 
following the classical fractal [16]. The semiclassical de- 
scription of probability density for such states has been 



derived recently in [17]. 

In view of the recent results described above I study nu- 
merically a simple model of the quantum Chirikov stan- 
dard map (kicked rotator) with absorption introduced 
in [18] which allows to vary continuously the fractal di- 
mension of the classical strange repeller. In this way 
the fractal Weyl law (1) is verified in the whole interval 
1 < d < 2. The model also allows to establish the lim- 
iting semiclassical distribution over escape rates 7 and 
find its links with the fractal properties of the classical 
strange repeller. The Chirikov standard map is a generic 
model of chaotic dynamics and it finds applications in 
various physical systems including magnetic mirror traps, 
accelerator beams, Rydberg atoms in a microwave field 
[19, 20, 21, 22]. The quantum model has been built up in 
experiments with cold atoms [23]. Thus the results ob- 
tained for this model should be generic and should find 
applications for various systems. 

The quantum dynamics of the model is described by 
the evolution matrix: 



i/> = Uip = Pe 



-iTh 2 /4 -ik cos 9 -iTn 2 /A 



(2) 



where ft = —id/dO and the operator P projects the wave 
function to the states in the interval [—N/2, N/2]. The 
semiclassical limit corresponds to k — > oo, T — > with 
the chaos parameter K = VP = const and absorption 
boundary a = N/k = const. Thus N is inversely pro- 
portional to the effective Planck constant T = h e ff, it 
gives the number of quantum eigenstates and the num- 
ber of quantum cells inside the classical phase space. The 
classical dynamics is described by the Chirikov standard 
map [19, 20] in its symmetric form: 



n + k sin 



Tn 
~2~ 



T 

1 + — (n + n) 



(3) 



Physically, the map describes a free particle propagation 
in presence of periodic kicks with period T (e.g. kicks of 
optical lattice in [23]). In this model all trajectories (and 
quantum probabilities) escaping the interval [—N/2, N/2] 
are absorbed and never return back. It is convenient to 
fix K = 7 so that the phase space have no visible stability 
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FIG. 1: (color online) Husimi functions of quantum fractal 
eigenstates with minimal value of 7 at N — 1025 (top left), 
N = 4097 (top right), N = 16349 (middle left) and the den- 
sity plot of classical strange chaotic repeller formed by or- 
bits nonescaping forward in time (middle right); two bottom 
panels show zoom for two middle panels respectively. Here 
a — N/k — 2, K — kT — 7 and the box counting dimen- 
sion of the repeller is d = 1.7230. In top and middle panels 
< 6 < 27r, -N/2 < n < N/2; density is proportional to 
color with red/gray for maximal density and blue/black for 
zero density. 

islands for a < 6 [24]. Then for the classical dynamics 
the probability P(t) to stay inside decays exponentially 
with time as P(t) ~ exp(— j c t) [18, 25], where j c is the 
classical escape rate and t is measured in the number 
of map iterations. For large values of a the spreading 
goes in a diffusive way and t c = l/j c ~ N 2 /D « 2a 2 
where D « k 2 /2 is the diffusion rate for K ^> 1. The 
independence of t c of N implies a = N/k = const. The 
quantum operator (2) can be considered as a simplified 
model of chaotic microlasers where all rays with orbital 
momenta below some critical value determined by the 
refraction index escape from a microcavity [9, 10, 11]. 

The right eigenstates and eigenvalues A m = 



exp(— ie m — 7 m /2) of the evolution operator U are de- 
termined numerically by direct dioganalization up to a 
maximal value N — 22001 (only states symmetric in n are 
considered). The Husimi distribution [26], obtained from 
smoothing of a Wigner function on a Planck constant 
scale, is shown in Fig. 1 for eigenstates with minimal 7 m 
at different values of TV at a = 2. With the increase of TV 
the Husimi distribution converges to a fractal set which 
is very similar to the classical strange repeller formed by 
classical orbits never escaping in the future. The classical 
repeller is obtained by iterating up to 3 x 10 9 classical tra- 
jectories homogeneously distributed in the whole phase 
space at t = 0. The classical remaining probability P(t) 
decays exponentially with 7 C = 0.2702 ±0.0011 and the 
computation of the box counting dimension [20, 21] of 
the strange repeller gives d = 1.7230 ± 0.0085. Accord- 
ing to [8, 20, 21] the information dimension d\ of the 
repeller can be expressed as d\ = 2 — 7 C /A, where A is 
the Lyapunov exponent. For large a and small j c it can 
be expressed via its value for the Hamiltonian dynamics 
on a torus where A « hi (if/ 2) = 1.2527 (for K = 7) [19]. 
This gives d\ = 1.7843 that is rather close to the nu- 
merical value of box counting dimension d (usually these 
two dimensions are rather close and, contrary to [13], I 
will not make difference between them). However, for 
smaller values of a the relation A = hi(K/2) is no more 
valid. To have A for all values of a its value is com- 
puted numerically following approximation used in [19]: 
A = (| ln(K cos(<9 + Tn/2)) |) where the average (...) 
is done over the orbits on the repeller. In this way A 
varies in the interval 1.913 < A < 1.294 for 0.7 < a < 6 
(A 1.363 and d 1 = 1.801 at a = 2). 

To check the validity of the fractal Weyl law (1) for 
various d the absorption border a is varied in the interval 
0.7 < a < 6 so that the classical fractal dimension and 
decay rate vary in the intervals 0.9976 ± 0.0060 < d < 
1.9367 ± 0.0067 and 1.6349 ± 0.0135 < 7 C < 0.0592 ± 
0.0003. In the quantum case the number of states N 7 is 
computed in the band width < 7 < 75 with 75 = 8/a 2 . 
In this way 75 > j c and the band contains a large fraction 
of fractal eigenstates. To improve the statistics, iV 7 is 
averaged over N r cases with slightly different values of fc± 
Sk with Sk < 2. Such a small variation of k does not affect 
the semiclassical properties but allows to improve the 
statistical accuracy. The number of realizations varied 
from N r = 40 at N = 101 to N r = 1 at N = 22001. 
The dependence of integrated number of states 7V 7 on 
N (x 1/h is shown in Fig. 2. The fit N 7 ex N v allows to 
determine the exponent v which according to (1) should 
satisfy the relation v = d — 1. It is important to stress 
that A^ 7 <C A^/2 at a < 2, so that the main part of 
(N + l)/2 eigenvalues has enormously large 7 ^> 75. 

The dependence of v on d is shown in Fig. 3. The 
law (1) is well satisfied for fractal dimensions 1 < d < 2. 
Certain deviations for d close to 1 should be attributed 
to rather small values of A^ 7 (e.g. N 7 = 26 at a = 0.7 and 
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FIG. 2: (color online) Dependence of the integrated number 
of states iV 7 with escape rates 7 < 7^ = 8/a 2 on matrix 
size A/". Symbols show numerical data for various values of 
absorption border a — N/k: 4 (circles), 2.5 (diamonds), 2 
(+), 1.5 (x), 1 (triangles), 0.8 (*), 0.7 (squares). Dashed lines 
show algebraic fits iV 7 oc N u with the fractal Weyl exponent 
v = 0.8930 ± 0.0028 (a = 4), v = 0.7129 ± 0.0073 (a = 2), 
1/ = 0.5697 ± 0.0042 (a = 1), v = 0.3559 ± 0.0094 (a = 0.8) 
v = 0.1175 ± 0.0069 (a = 0.7, here 7^ = 4/a 2 ). Logarithms 
are natural. 



N = 22001) so that even larger N values are required to 
see the asymptotic behavior. The relation v = 1 — 7 C /A 
also works rather well even if the Lyapunov exponent A 
should be probably computed in a more exact way for 
small values of a ~ 0.7. Thus, the data of Fig. 3 confirms 
the validity of the fractal Weyl law in the whole available 
interval of fractal dimensions. 

In addition to the integrated characteristic (1) it is in- 
teresting to consider the differential distribution dW/dj 
which determines the number of states in the interval dj 
at given 7. The evolution of distributions dW/dj with 
the growth of N is shown in Fig. 4. The data clearly 
show that in the semiclassical limit dW/dj converges to 
a certain limiting distribution independent of N. Such an 
effect has been noticed already in first studies [18] where 
mainly the diffusive limit with a = 10 and ) c « A has 
been considered. In such a case the dimension is very 
close to the integer value d — 2 and due to that the frac- 
tal dependence (1) has been missed in [16, 18] even if the 
fractal structure of eigenstates has been clearly detected 
[16]. In the diffusive case d « 2 one has dW/dj oc I/7 3 / 2 
for 7 > 7c that is explained by simple estimates [18] and 
more rigorous analytical treatment [27]. When the frac- 
tal dimension d is noticeably less than 2 than j c ~ A and 
the diffusive approximation is no more valid. A distinc- 
tive feature of the distribution in this case is the gap in 
the distribution dW/dj which is zero for 7 < 7 C , sharp 
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FIG. 3: (color online) Dependence of the fractal Weyl expo- 
nent v on the fractal box counting dimension d: full circles 
show numerical data, the straight line shows the fractal Weyl 
law (1) v — d — 1, pluses show 1 — 7 C /A versus d — 1 which 
should follow the relation di — 1 = 1 — 7c/A, where A is the 
Lyaponov exponent computed approximately (see text). Here 
0.7 < a < 6, 0.1175 < v < 0.9402, 0.9976 < d < 1.9367. 
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FIG. 4: (color online) Dependence of the distribution 
jcdW/dj on the rescaled escape rate 7/7C for different values 
of N at a = 2. Here 7 C = 0.2702, d = 1.7230, the probability 
dW/dj is normalized to unity in the interval jb and N is 22001 
(blue/black full curve iV 7 = 1278), 12801 (maroon/gray full 
iV 7 = 1022), 6401 (green/gray full iV 7 = 500), 3201 (red/gray 
full Nj = 293), 1601 (orange/gray dashed iV 7 = 181), 801 
(turqse/gray dashed iV 7 = 114), 401 (magenta/gray dashed 
iV 7 = 68.7). 



peak at 7 = 7 C followed by a smooth drop at 7 > j c (this 
drop is compatible with dependence I/7 3 / 2 ). 

These properties of the distribution dW j dj remain es- 
sentially the same when j c is changed by a factor 3.5 as it 
is shown in Fig. 5. Indeed, the shape of the distribution 
varies very little for 1.5 < a < 4 and becomes broader 
only at a < 1.5 . The later case have however relatively 
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FIG. 5: (color online) Same distribution as in Fig. 4 drown 
for various absorption boundaries a at N = 22001. Here a is 
4 ( 7c = 0.1019, iV 7 = 3607, black curve); 2.5 ( 7c = 0.2063, 
iV 7 = 2032, orange/gray curve); 2 ( 7c = 0.2702, N 7 = 1278, 
blue/black curve); 1.5 ( 7c = 0.2961, iV 7 = 1342, green/gray 
curve); 1 ( 7c = 0.6967, iV 7 = 472, red/gray dashed curve). 

small statistics iV 7 and probably larger N should be used 
to reach a limiting distribution for a < 1.5. It is interest- 
ing to note that dW/dj has certain similarities with the 
Wigner proper times distribution discussed in [28]. 

In conclusion, the obtained data confirm the validity 
of the fractal Weyl for all fractal dimensions in the in- 
terval 1 < d < 2. They show the existence of the limit- 
ing distribution of the Gamow resonances dW/dj which 
has a gap of size 7c above which the distribution has a 
sharp peak (see Figs. 4,5). Thus the classical decay rate 
7c essentially determines the quantum decay rates on the 
quantum fractal corresponding to the classical strange re- 
peller with orbits never escaping in future times (Fig. 1). 
The analytical computation of the limiting distribution 
dW/ dj still remains an open problem. It is possible that 
the analytical methods pushed forward recently [17] will 
allow to make progress in this direction. Also, it would 
be interesting to check the validity of the fractal Weyl 
law for dimensions d > 2. In such a case it is natural to 
expect that 7V 7 ex Ti~ v with v = d — rif where d is the 
fractal dimension of the classical strange repeller and rif 
is the number of degrees of freedom (in the present model 
rif = l,d< 2). 

At present the properties of large nonunitary matri- 
ces find important applications in various areas includ- 
ing search on the Internet [29, 30] and it is possible that 
the fractal quantum eigenstates may have there certain 
applications since they give an example of important non- 
trivially connected fractal sets of small measure. 

I thank S. Nonnenmacher for highlighting of his results 
and K.M. Frahm for stimulating discussions. 
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